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Abstract For the nonthermal radio emission of the Galactic Center Arc in situ electron acceleration is imperative. 
The observed radio spectrum can be modeled by a transport equation for the relativistic electrons which includes 
particle acceleration by electric fields, momentum diffusion via scattering by magnetohydrodynamical turbulence 
and energy losses by synchrotron radiation. The accelerating electric fields can be regarded as a natural conse- 
quence of multiple reconnection events, caused by the interaction between a molecular cloud and the Arc region. 
The radio spectrum and even the recently detected 150 GHz emission, explicitely originating from the interac- 
tion regions of a molecular cloud with the magnetized Arc, can be explained in terms of quasi-monoenergetically 
distributed relativistic electrons with a typical energy of about 10 GeV accelerated in stochastically distributed 
magnetic reconnection zones. 
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1. Introduction 

The Galactic Center (GC) Arc is a unique nonther- 
mal radio continuum structure in the Galaxy. It is lo- 
cated at a projected distance of about 30 pc from 
the Sgr A complex and consists of a network of mag- 
netic filaments filled with relativistic electrons running 
almost perfectly perpendicular to the Galactic plane 
(e-g- 



(Yusef-Zadeh et al. 1984) I. The nonthermal nature 



of these filaments has been proven beyond any doubt by 
radio emission which exhibits a degree of linear polariza- 
tion which is close to the intrinsic value at high frequen- 
cies of 60% |(Reich et al. 1988)|[(Lesch fc Reich 1992)| We 
emphasize the nonthermal character of the Arc since its 
radio spectrum is a very remarkable one. A decomposi- 
tion of its radio spectrum between 843 MHz and 43 GHz 
by (Reich et al. 1988) shows a spectral index a = +0.3 
(we use the conventional relation that the observed flux 
scales as i^"). This finding of an increasing radio spec- 
tral index has been confirmed by high resolution VLA- 



observations (Anantharamaiah et al. 1991) 



Such an inverted radiation spectrum with an index 
of -1-0.3 is expected from a quasi-monoenergetic elec- 
tron distribution or an energy distribution with a well- 



defined low energy cutoff, respectively (Lesch et al. 1988) 
Observations at high frequencies (between 32 GHz 
(Lesch fc Reich 1992) | and 43 GHz |(Sofue et al. 1987b) I 
seem to indicate a spectral turn over, i.e. a fading of 
the Arc towards higher frequencies. However, Reich et al. 
(2000) detected an enhanced emission at 150 GHz slightly 



offset relative to the most intense vertical nonthermal fil- 
aments seen at lower frequencies. This emission originates 
from the apparent interacting areas of dense molecular 
material with the Arc. 150 GHz emission was also detected 
south of a molecular cloud, where the vertical filaments of 
the Arc cross a weak filamentary structure. The spectrum 
of the emission is inverted relative to 43 GHz and compat- 
ible with an origin from quasi-monoenergetic electrons or 
an electron distribution with a low energy cutoff, but not 
with optically thin emission from cold dust. Reich et al. 
(2000) conclude "that the coincidence of enhanced emis- 
sion with regions of interacting molecular gas strongly sug- 
gests that high-energy electrons are accelerated in those 
places where the magnetic field is compressed" . They cal- 
culated a Lorentz factor of the electrons to be 7 ~ 2 x 10'' 
emitting synchrotron radiation at 150 GHz within a mag- 
netic field of ImG. 

To summarize, in the Galactic center several areas ap- 
pear to be filled with relativistic electrons whose energy 
must be distributed quasi-monoenergetically around a few 
GeV. We note that even the very center of our Galaxy Sgr 
A* exhibits a radio to infrared spectrum which is in sur- 
prisingly good agreement with optically thin synchrotron 
emission of a quasi-monoenergetic electron distribution 
with a typical energy of 80 MeV, i.e. Lorentz factor of 
160 



(Duschl fc Lesch 1994; 



In a former paper about particle acceleration in the 
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Arc (Lesch & Reich 1992) we considered the magnetohy- 
drodynamical interaction of a molecular cloud with the 
Arc filaments. Especially we investigated the role of the 
moving gas cloud as a trigger mechanism for magnetic field 
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amplification accompanied by dissipation of the magnetic 
energy via magnetic reconnection. 

Magnetic reconnection takes place when magnetic field 
lines with antiparallel directions encounter. Such a situa- 
tion corresponds to the formation of an electric current 
sheet. Typically astrophysical plasmas are ideal electric 
conductors, i.e. their electrical conductivity is very high. 
In such media the magnetic field is frozen into the motion 
of the conducting fluid. Any plasma velocity distortion 
onto the magnetic field lines is automatically related to an 
electric field which is perpendicular to the plasma velocity 
and the magnetic field. This electric field is described by 
the ideal Ohm's law 



f 



vx B = 0. 



(1) 



In a plasma where the magnetic field is strongly distorted 
by shear flows, radial explosive flows or stochastic mo- 
tions, it is unavoidable that fleld lines with antiparallel 
directions encounter. Consequently, the ideal conducting 
plasma switches into a nonideal medium with a flnite, lo- 
calized electrical conductivity, i.e. the ideal form of Ohm's 
law is violated. The strong spatial gradients of the mag- 
netic fleld in such interaction zones represent an energetic 
" crisis" which is relaxed by partial dissipation of the mag- 
netic energy via the formation of current sheets in the 
resistive medium (e.g. (Priest & Forbes 2000) and refer- 
ences therein). The nonideal Ohm's law is given by 



1 



-V X B = R^O, 



(2) 



where R is some yet unspecifled nonideal term, i.e. the 
plasma resistance. 

In other words, magnetic reconnection corresponds 
to a magnetic fleld aligned electric fleld which 
is associated with a generalized electric potential 
(Schindler et al. 1988)| [(Schindler et al. 1991)| 



V 



Es ds = 



Rs ds 



(3) 



integral is evaluated along the mag- 
lines that penetrate the reconnection 
such a potential drop offers 



where the 
netic fleld 
region. Of 
the possibility to accelerate 
ciently {( Schindle r et al. 1991) 



course, 



(Lesch fc Birk 1997)t 
(Schopper et al. 1998)[ 



particles very effi- 
KBlackma n 1996)! 



(Lesch fc Birk 1998)1 



(Nodes et al. 2003) 



(Litvinenko 1999)| 



Instead of investigating the acceleration in one re- 
connection zone, we consider in this contribution the 
acceleration of relativistic electrons in numerous recon- 
nection regions, i.e. current carrying fllaments driven by 
the interaction of a molecular cloud with the magnetic 



fleld in the Galactic Center Arc (Lesch fc Reich 1992) 



(Serabyn & Morris 1994) This scenario arises quite natu- 
rally, since the necessary energy source is represented by 
moving molecular gas which encounters the poloidal mag- 
netic flelds in the Arc which is proven to be present by 



the observed very high polarization of the radio emission 
up to 60%. 

Since we have many acceleration regions we apply a 
dynamical description of the energy distribution function 
of relativistic electrons based on systematic momentum 
gains by reconnection, losses by synchrotron losses and 
momentum diffusion due to scattering on magnetohydro- 
dynamical turbulence, i.e. Alfven waves and magnetosonic 
waves. 



2. The temporal evolution of the energy 
distribution function of relativistic particles 

The dynamical evolution of a population of high-energy 
particles that do not interact with each other is described 
by the Liouville equation 



dF dr„ dp„ 



(4) 



where F(r, p, t) is the one-body distribution function. Due 
to the complexity of the non-linear dynamics in 6-1-1 di- 
mensions this equation is only of limited use. Fortunately, 
for our purposes we do not need to know the full informa- 
tion of the particle dynamics. Rather, we will make use 
of the following assumptions. First, we will assume pitch 
angle isotropy. This should be granted by efficient parti- 
cle scattering in momentum space in reconnection regions 
caused by Alfvenic and magnetosonic wave flelds with 



relatively high-energy densities (see (Schlickeiser 1986) I. 
Second, we are interested in the global energy spectrum 
of the radiating electrons and thus, can dispense with the 
detailed spatial dependence of the distribution function. 

Consequently, we deal with the isotropic distribution 
f{p,t) — J drF{r,p,t)/V where V denotes the emission 
volume. The total particle number is calculated as N{t) = 
An J fp^dp. The continuous momentum gains and losses 
Ap of the particles can be described by the momentum 
operator 



9 [l,Pgain I Ploss 

p"^ op 



)p' 



(5) 



where the systematic temporal momentum gains and 
losses are denoted by Pgain and pioss, respectively. 

Additionally, we will allow for stochastic diffusion in 
momentum space by scattering of magnetohydrodynam- 
ical fluctuations (in the sense of a Fokker-Planck term) 
denoted by Dp. In principle, catastrophic losses that may 
result to a further sink in the balance equation for / can 
be modeled by a characteristic time r, i.e. df /dt ^ —J/t. 
Finally, the injected particle population is represented by 
some explicit source term S{p). Accordingly, the tempo- 
ral evolution of the particle distribution is governed by the 



equation,(Schlickeiser 1984) (Schlickeiser 1986) 



dfjp^t) 1 d / ^^ dfjp^t) 
dt p^dp\ ^ dp 



■'Cpf{p,t)- 



f{p,t) 



s{pm 
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which holds, if the distribution function f(r,p,t) is sep- 
arable in space and momentum, and if the momentum 
operator Cp does not depend on the spatial coordinates. 
In the present application based on the reconnection sce- 



detected cloud motion is in the range of 15-45 km/s 



nario (see also (Birk et al. 2001) I the momentum gain is 
given by Pgain = diil'^^) = l^s = C = const., where 7 
denotes the Lorentz factor, Eg is the accelerating electric 
field component (locally directed parallel to the magnetic 
field) and q is the electrical charge. 

We emphasize that the acceleration of high-energy par- 
ticles is regarded as a consequence of the conversion of 
magnetic field energy into particle energy by magnetic re- 
connection. Since synchrotron losses are responsible for 
the observed non-thermal radio emission from the Arc 
we take pioes = —kp"^ with k — ^ctt^Ub / m^c^ , with the 
Thomson cross section ctt 1 the electron mass rric, the 
speed of light c and the magnetic field energy density Uq . 
For our applications inverse Compton scattering, which 
has the same 7-dependence as synchrotron losses, as well 
as bremsstrahlung are negligible. In case of particle scat- 
tering in momentum space by Alfvenic waves the momen- 
tum diffusion coefficient is given (Schlickeiser 1984) by 
Dp = Dp^ with D = v\/9K{p), where va is the Alfven 
velocity and K{p) — Xc/3/3 is the spatial diffusion coeffi- 
cient, which has been expressed in terms of the particle's 
mean free path A [(Sch lickeiser 1986) For our purpose we 
have set f3 = v/c = 1. Since we have to model the forma- 
tion of monoenergetic distribution functions, the effect of 
catastrophic losses must be negligible, because otherwise 
the particles would leave the acceleration region before a 
pile up in energy has been established , (Schlickeiser 1986)| 

3. Numerical results 

In order to calculate the temporal evolution of f{p,t) by 
numerical integration of equation 6 we first have to spec- 
ify the physical parameters k, D, and S(j)), which are 
related to the physical conditions of the Galactic Center 
Arc. First of all, we discuss the source term S{p). Since 
the energy distribution of the injected particles, which are 
generated in the Galactic center, is a monoenergetic one 
(Lesch & Reich 1992) the use of a (5-like source term, rep- 
resented, e.g., by a cos/i~^ (p)-function in our calculations, 
is a reasonable assumption. As a typical injection energy 
we use 75 MeV, the value close to the one calculated by 
(Duschl fc Lesch 1994)|for the central object Sgr A*, i.e. 



7„j = 150. 

The magnetic field strength in the arcs fila- 
ments is of the order of a few mG, or higher 



(Yusef-Zadeh & Morris 1987a) To be on the save side we 
use the 150 GHz detection of (Reich et al. 2000) close to 
the Arc as the maximum synchrotron frequency emitted. 
Thus, we have to deal with a Lorentz factor of about 
7 ~ 2 • 10'', necessary for 150 GHz synchrotron emission 
in a magnetic field of 10~^ G. 

The maximum electric field strength is given by 
E ~ w/cB, where w denotes the bulk velocity of the 
thermal plasma that acts as an MHD generator. The 



(Tsuboi et al. 1997) Again to be on the save side we will 
use w = 5 km/s which gives a maximum electric field 
strength of 10~^ statvolt cm^^ . If we assume, that elec- 
trons can be accelerated over the whole extension of the 
Arc (about 40 pc) without energy losses (ideal linear ac- 
celerator), we need a minimum strength of the average 
electric field of i?s — 3 • lO^^'^ statvoltcm^^ to reach the 
observed particle energies, which is a tiny fraction of the 
available electric field. 

This oversimplification is helpful to work out the anal- 
ysis with a reference value of C and is not meant as a 
physical assumption we apply (see discussion below). 

If we restrict our calculations to electrons we obtain 
an estimate for the systematic acceleration term 



Pgain = e^^s - C > 1-5 • 10-22 dyne 



(7) 



Due to the high polarization of the synchrotron radiation, 
up to 60 % (Lesch & Reich 1992) we consider only syn- 
chrotron emission to be responsible for the energy losses 
of the particles in the Arc. In this case we can find a 
certain value po = [SCm'^c'^ /AaTllBy^^ where momentum 
gains and losses keep the balance. For momenta p > po 
momentum loss dominates momentum gain. For momenta 
p < Po the opposite holds. Thus, we expect a pile up for 
the distribution function f(j),t) at po, which is consistent 
with a monoenergetic d istribution of the electrons (see also 
(Lesch & Reich 1992) I. As discussed above, the pile up 
in the Arc should occur at momenta with Lorentz factor 
7 ~ 2 • 10**. If one assumes, for simplicity, a constant elec- 
tric field of -Eg ~ 3 • 10-^3 statvolt cm ^ over the complete 
expansion of the arc, a magnetic field of 3 • 10"'^ Gauss is 
necessary for a pile up at po = JtUcC — 5.4 • 10^^'^ g cms"^ 
according to a systematic loss coefficient of 



k — — ctt 



1.4 • 10^ g" 



(8) 



However, the particles, in fact, will be accelerated locally 
in numerous reconnection regions with £'|| 7^ and will 
be scattered by magnetohydrodynamic fluctuations like 
Alfvenic and/or magnetosonic waves. Consequently, the 
electric field strength in the reconnection regions must be 
ultimately stronger than 10~^^ statvolt cm"-'. Indeed, the 
electric field can be estimated by the ratio of the number 
density of relativistic electrons rir and the thermal electron 



density (Papadopoulos 1977) 



exp 



E^ 
2E 



(9) 



Ec — TOe^thet'coii/e dcnotcs the critical electric field, 
whereas z;the is the thermal velocity of the electrons and 
t'coii is the collision frequency. As calculated by Lesch & 
Reich 1992 the thermal electron number density ric in the 
Arc filaments is about 10 cm^"^ and the density of rela- 
tivistic electrons is of the order of 10^^ cm^^^. If we use an 
electron temperature of 10^ K we obtain from equation 9 
E o± 3 ■ 10^^° statvolt cm~^, which can be easily achieved 
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by the available electric field, induced by the w x B - 
motion of a molecular clouds calculated above. 

To obtain the diffusion coefficient D we first have to 
estimate the particle's mean free path A. As shown in 
(Schlickeiser 2002)1 A can be determined by 



where vq is the gyro radius of the protons, Bq/SB gives 
the ratio between macroscopic and fluctuated magnetic 
field strength, Lq is the expansion of the considered ob- 
ject and q is turbulence spectral index. In the vicin- 
ity of a strong magnetic distortion, like the Arc fil- 
aments. Burgers turbulence is expected to be excited 



(Chambers et al. 1988) Thus, we choose a turbulence 
spectral index of g = 2. With a particle density of n ^ 
10 cm~^ and a magnetic field strength of 3 • 10"'^ Gauss 
the Alfvenic velocity is about 700 km s~^. Hence we ob- 
tain for the diffusion coefficient 



D = 



3cA 



SB^ 

Bo 



2 2 
cLq 



5.4-10 



-18„ 



(11) 



if we assume, that Bq/SB is of the order of 10 which im- 
plies a reasonable level of well-developed turbulence. Much 
higher diffusion coefficients would inhibit pile-up distribu- 
tions. However, a quasi-linear theory applied for equation 
9 requires SB < Bq. 

The synchrotron emission spectrum associated 
with the particle spectrum can be calculated from 



(Rybicki and Lightman 1979) 



where 



PMN{p) dp. 



~3 



^|(0d?, 



(12) 



(13) 



is the power per unit frequency v emitted by each electron. 

5 (^) denotes the modified Bessel function of the second 
kind and = 'A^'^loq/ At: is the critical frequency {ujq is 
the electron gyro frequency). 

Having clarified the initial conditions and chosen phys- 
ical parameters we now come to the numerical results. 
Figure 1 shows the temporal evolution of the particle 
density iV(7) and the corresponding synchrotron emission 
spectrum, calculated from equation 12. In the chosen nor- 
malization the time is measured by Iq = 10^ s. As long as 
the particle acceleration exceeds the synchrotron losses, 
the quasi-monoenergetic distribution of the electrons is 
shifted toward higher momenta (a, c). Due to momen- 
tum diffusion the maximum of N{'^) decreases with time 
(c). As soon as the synchrotron losses are comparable to 
the energy gains, the maximum of N{'^) increases again 
(Figure 2 (a)). At a Lorentz factor of about 7 = 2- 10^ the 
momentum gains and losses are balanced and the distri- 
bution of the electrons that are assumed to be contained 



in the Arc during the entire radiation processes devel- 
ops towards a shifted pile up. Consequently, a monoen- 
ergetic distribution function evolves again (Figure 2 (c)). 
The cut-off frequency remains almost constant during en- 
tire pile-up dynamics. The associated radiation spectrum 
(Figure 2 (c)) corresponds to the one observed by sev- 
eral authors, including the 150 GHz emission detected by 



(Reich et al. 2000) We note that for a finite loss time of 
T < 10^*^ s the pile-up is prevented by the particle losses. 

4. Conclusions 

In the filaments of the Galactic Center Arc electrons are 
accelerated to considerably high energies of about 10 GeV. 
This is obvious from the detected polarized radio emission 
exhibiting a very high degree of polarization (60%) and a 
rising spectrum with a spectral index of +0.3 up to 150 
GHz. Such a spectral behavior is due to an energy distribu- 
tion function of the radiating relativistic electrons which 
either has a low-energy cutoff or which is quasi mononen- 
ergetic. No obvious energy sources for particle accelera- 
tion are present in the neighborhood of the filaments. The 
very center of the galaxy, SgrA* is a source for monoener- 
getic relativistic particles but on significant lower energies 



of about 50 MeV |(Duschl fc Lesch 1994)| However, the 
magnetic filaments interact with the plasma of a molecu- 
lar cloud which moves with velocities of some 10 km s~^ 
relative to the magnetic field. This interaction can be in- 
terpreted in terms of the induced Lorentz force inducing a 
convective electric field which in the first place is oriented 
perpendicular to the magnetic field and the local plasma 
velocity of a molecular cloud. The motion of the plasma lo- 
cally distorts the magnetic field necessarily in such a way 
that antiparallel directed field lines encounter. Provided 
that some violation of ideal Ohm's law occurs, e.g. due to 
microturbulence, magnetic reconnection will occur. The 
magnetic field energy will be partly converted to parti- 
cle energization. In the reconnection regions a magnetic 
field-aligned electric field component forms with a magni- 
tude of some fraction of the convective electric field. This 
parallel electric field represents a perfect candidate for ef- 
ficient particle acceleration. The required electric field for 
the electron energization are considerably weaker than the 
magnetic field which is in accordance with the observa- 
tional fact that the cloud plasma moves with velocities 
much lower than the speed of light. Since reconnection 
is a localized phenomenon which depends on the local 
properties of the plasma, we consider a scenario charac- 
terized by multiple reconnection sheets in which particles 
are accelerated and scattered by magnetohydrodynamical 
fluctuations. Such scattering leads to diffusion in energy 
space. Moreover, we take into account the energy losses by 
the observed synchrotron radiation. By means of a rela- 
tively simply transport equation for the distribution func- 
tion of the radiating relativistic electrons we could show 
that an injected monoenergetic distribution function with 
75 MeV (coming from the Galactic center) evolves into 
a quasi-monoenergetic distribution function at 10 GeV. 
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The resulting radiation spectra increase with frequency 
up to some hundred GHz according to the observations, 
i.e. ^ 

A major assumption of our study is the complete 
isotropy of the distribution function which is based on 
calculations by (Achatz et al. 1990) who show that the 
isotropisation time tiso = c/vAiB/dB)l/n^ where fl is the 
gyro frequency of the protons. In this application ^iso is 
about 4 • 10^ s which is much shorter than the studied 
time scale of the evolution of the distribution function. 
Some non-isotropic particle population would lead to an 
overlap of a second radiation component in a specific en- 
ergy range that we can not handle in our model. However, 
the present observations seem to give no hint to such a 
component. 
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(a) Temporal evolution of 
A'' (7). The solid, dotted and 
dashed lines show the injected 
particle spectrum and the evo- 
lution after t = Ito and t = 
2 to, respectively. 
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(b) The corresponding syn- 
chrotron spectrum I^, 
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(c) Temporal evolution of 
A'' (7). The solid and dotted 
lines show the evolution af- 
ter t = 10 to and t = 20 to, 
respectively. 



(d) The corresponding syn- 
chrotron spectrum 7„ 



Figure 1. Temporal evolution of A''(7) calculated from 
equation 6. The values of the physical parameters, which 
we have used for our calculations are: ( = 1.5 • 10~^^ dyne, 
k = 1.4-103 g-icm~i and D = 5.4-10-^Ss-^ To evaluate 
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(a) Temporal evolution of 
N{-y). The sohd and dotted 
lines show the evolution after 
t = SOto and t = 50to, respec- 
tively. 




(b) The corresponding syn- 
chrotron spectrum 
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(c) Temporal evolution of 
A/'(7). The solid, dotted and 
dashed lines show the evolu- 
tion after t = 80to, t = 100 to 
and t = 150 to, respectively. 




(d) The corresponding syn- 
chrotron spectrum 7„ 

Figure 2. Temporal evolution of iV(7) towards a quasi- 
monoenergetic distribution function (a), (c). The result- 
ing synchrotron emission spectrum (d) shows an inverted 
spectrum with a spectral index of a = 1/3 {Ii, ~ v"") and 
a cutoff-frequency of some hundred GHz after 1.5 • 10^^ s. 



